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ABSTRACT 

Stability properties of "magnetic tower" jets propagating in the gravitationally stratified background have 
been examined by performing three-dimensional magnetohydrodynamic simulations. The current-carrying, 
Poynting flux-dominated magnetic tower jet, which possesses a highly wound helical magnetic field, is subject 
to the current-driven instability (CDI). We find that, under general physical conditions including small pertur- 
bations in the initial background profiles, the propagating magnetic tower jets develop the non-axisymmetric, 
m = 1 kink mode of the CDI. The kink mode grows on the local Alfven crossing time scale. In addition, two 
types of kink modes appear in the system. At the central region where external thermal pressure confinement 
is strong, only the internal kink mode is excited and will grow. A large distance away from the central region 
where the external thermal pressure becomes low, the external kink mode is observed. As a result, the exterior 
of magnetic tower jets will be deformed into a large-scale wiggled structure. We also discuss extensively the 
different physical processes that contribute to the overall stability properties of the magnetic tower jets. Specif- 
ically, when the jet propagates in an initially unperturbed background, we find that they can survive the kink 
mode beyond the point predicted by the well-known Kruskal-Shafranov (K-S) criterion. The stabilization in 
this case comes mainly from the dynamical relaxation of magnetic twists during the propagation of magnetic 
towers; the magnetic pitch is reduced and the corresponding K-S critical wavelength becomes longer as the 
tower jet proceeds. Furthermore, we show that the pressure-driven and Kelvin-Helmholtz instabilities do not 
occur in the magnetic tower jets. This strongly suggests that the CDI is the primary reason for the wiggling 
structures in jets. 

Subject headings: instabilities — galaxies: active — galaxies: jets — methods: numerical — MHD 



1. INTRODUCTION 

Magnetohydrodynamic (MHD) mechanisms are often in- 
voked to explain the launching, acceleration and collimation 
of jets from Young Stellar Objects, X-ray binaries, Active 
Galactic Nuclei ( AGNs ), Microquasars, and Quasars (see, 
e. e.. iMeier et al. 1 1200 It and references therein). Strongly 
magnetized jets, particularly those with a strong toroidal 
field encircling the collimated flow, are often referred to as 
"current-carrying" or "Poynting flux-dominated" (PFD) jets. 
A large current flowing parallel to the jet flow is responsible 
for generating a strong, tightly wound helical magnetic field. 
The global picture of a current-carrying jet with a closed cur- 
rent system linking the magnetosphe re of the central engine 
and the hot spots was introduced by Benford ( 1978. 120061) 
and applied to AGN double radio sources. This closed cur- 
rent system includes a pair of current circuits, each containing 
both a forward electric current path (the jet flow itself, with its 
toroidal magnetic field, toward the lobe), and a return electric 
current path (along some path back to the AGN core). 

Recent high resolution observations, in X-ray, opti- 
cal and radio bands, show that many AGN and quasar 
jets display wiggles/kinks, which may be associated 
with a helical motion of the underlying bulk flow, in 
wide spatial ran g es from sub-parsec t o kpc scales (e.g. , 
Reid et al. I fl989t iHummel et al. I fl99l iBaum et all fl997l 
Feretti et al. ll999ULobanoy & Zensusl2001tlMarshall et al. I 
20011 lEilek & Owen 1120021 iSavolainen et al. II2006I) . Many 
ideas have been proposed for these helical distortions: MHD 
instabilities such as the Kelvin-Helmholtz instability (KHI), 
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the current-driven instability (CDI), and the pressure-driven 
instability (PDI), or the precession of the jet ejection axis 
due to the existence of a binary black hole ( Beg elman et al. I 
1980), or an encounter with another galactic core. 

It is well known that a cylindrical plasma column with a he- 
lical magnetic configuration is subject to the m = (sausage) 
mode of PDI, the m = 1 (kink) mode of CDI, and the other 
higher order modes, where m is the azimuthal (4>) mode num- 
ber (s ee, e. e. . iKadomtsev ifl 966t iBateman |[l978t iFreidberg I 
119821 for details). CDI is driven by parallel (axial) current 
flows and can exist even in a zero pressure, force-free plasma. 
On the other hand, PDI is powered by perpendicular cur- 
rent flows and it is often known as the interchange (internal) 
modes. KHI is considered to be very important at the shearing 
boundary between the flowing jet and the external medium, 
particularly in the kinetic flux-dominated (KFD) jets. In prin- 
ciple, CDI and KHI have both external (surface) and internal 
(body) modes. PFD jets should be especially susceptible to 
CDI rather than PDI and KHI. This is because the strong ax- 
ial electric current is present inside the PFD jets, which is 
responsible for the highly wound helical field configurations, 
and the Lorentz force plays an important role in the jet dy- 
namics. The investigation of the destructive influence of CDI 
on PFD jet flows is an important avenue of research in the 
study of astrophysical jets. 

Much more attention on studies of stabilities (disruption of 
the interior/exterior of astrophysical jets) have given to KHI 
than to CDI. This is because jets were believed to be super- 
Alfvenic or super-fast magnetosonic KFD flows {i.e., kinetic 
energy was expected to greatly exceed the magnetic energy). 
However, it is not necessarily always the case; a super- 
Alfvenic or super-fast magnetosonic PFD fl ow can exist in 
both n onrelativistic (Kudoh & Shibata 1997; Vlahaki s" et al. I 
2000) and relativistic lviahakis & Konigl 11200 ll) MHD out- 
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flow solutions. Analytic studies on the CDI as a possible 
explanati on for Jet disruption have been p erformed as well 
( Eichl erl[T993l ISpruit. Fodizzo. & Stehle 1119971 iBegelman 
1998; Lvubarskii 1999). Linear analyses of nonrelativis- 
tic force-free jets, which are thermally confined by the ex- 
ternal n onmagn etized medi um, were per formed in sup er-fast 
JAppI & Camenzind1ll992l) or trans-fast (lAppll [l996) mag- 
netosonic regimes. They found that the growth rate of CDI 
is always substantially smaller that those of KHI in their spe- 
cific treatment. Only for quite small fast magnetosonic num- 
bers, Mp < 0.5, both become comparable. Linear anal- 
yses of the relativistic force-free jets, in case of constant- 
longitudinal field B? along the jet axis, were performed by 
Istomin & Parievl (Tl 9941 1 1 9961) . They concl uded that such a 
jet is stable against CDI, while Lvubarskii (1999) analyzed 
the relativistic force-free jets with B z decreasing outwards 
and found that they are unstable. 

Nonlinear development of the C DI for cold super-fa st mag- 
netosonic jets has been studied by Lerv et al. (2000), based 
on the linear analysis (Applet al. 2000]). It was found that 
the current density is redistributed within the inner part of the 
jet radius on a characteristic time scale of the order of the 
Alfven crossing time in the jet frame. Nothing in their nu- 
merical results indicated a possible disruption of the jet by 
the CDI sausage (m = 0) or kink (m = 1) modes (see also 
iLerv & Frank l2000tlFrank et al. 120 00). In general, these lin- 
ear/nonlinear considerations conclude that "current-driven in- 
stabilities are unli kely to be da ngerous for the integrity of as- 
trophysical jets" (lAppll fl996) and this has been widely ac- 
cepted in theoretical/observational communities. In the non- 
linear regime, an interplay between KHI and CDI can occur 
in the helically magnetized flo w even with slightly sup er-fast 
magnetosonic speed Mp > 1 (Batv & Keppens 2002). This 
nonlinear interaction can contribute to jet survival, and the 
large-scale magnetic deformations associated with the CDI 
mode development can effectively saturate KHI surface vor- 
tices and prevent jet disruption. 

By contrast, the first attempt to investigate the nonlinear 
behavior of visible distortion in KFD jets along the large- 
sc ale magnetic field by the CDI kink mode was performed 
bv iTodo et al. I (119931) to mode l Herbig Haro objects. In a 
similar study, Nakamura et al. (2001) examined the forma- 
tion of wiggled structures by the m = 1 mode in PFD jets to 
apply for AGN jets. Their numerical result is applicable for 
the system atic distribution of the observed rotation measure 
in 3C 4 49 (iFeretti et al. |[T9 99). associated with the jet defor- 
mation (Kigure et al. 1 120041) . The rotating PFD/KFD jets can 
be potentially stabilized against the m = 1 mode beyond the 
original Kruskal-Shafranov criterion, however, sudden desta- 
bilization of rotating jets due to the rn = 1 mode can occur via 
the angular momentum lo ss triggered by MHD shock waves 
JNakamura & Meier 12004 . 

The purpose of the present paper is to discuss the sta- 
bility properties of "magnetic tower" jets, which have not 
been studied in the past work. Of particular interest here 
is the nonlinear growth of the CDI kink mode. The CDI 
kink mode in magnetic tower jets may be distinguished 
from that in the classical MHD j e ts along the large-scale 
magn etic field (|Todo et alH 119931 iNa kamura et al. 2001; 
Nakamura & Meier 2004), while it should be co mparable di- 
rectly to_||n2n^srupdng" models (Appl 1996; App l et alTI 
1200(1 ILerv et al. II2000I) . in which the MHD jets are assumed 
to be in pressure equilibrium with an unmagnetized ambient 
medium. 



This is the third of our series of papers to examine the non- 
linear magnetic tower jets. The first of our series. IL1 et al. I 
(2006, hereafter Paper I), described the basic assumptions 
and the approaches in the numerical modeling of magnetic 
tower jets. The evolution of the tower jets in a constant den- 
sity/pressure bac kground has been exam ined there. The sec- 
ond of our series, Nakamura et al. ] (120061 hereafter Paper II) 
investigated the global structure of magnetic tower jets in the 
gravitationally stratified atmosphere (a more realistic astro- 
physical situation), in terms of the MHD waves structure, the 
radial force equilibrium, and the collimation. In the present 
paper, the third of our series, we investigate the stability prop- 
erties of magnetic tower jets. We examine our numerical re- 
sults in two cases; one without initial perturbations and an- 
other with initial perturbations. The former exhibits the quasi- 
axisymmetric evolution without any growth of MHD instabil- 
ities, while the latter shows non-axisymmetric evolution by 
the CDI modes. This paper is organized as follows. In §2, 
we outline our numerical methods. In §3, we describe our 
numerical results. Discussion and conclusion are given in §4 
and §5. 

2. NUMERICAL METHODS AND MODEL ASSUMPTIONS 

We solve the nonlinear system of time-dependent ideal 
MHD equations numerically in a 3-D Cartesian coordinate 
system (x, y, z). The basic numerical treatments (including 
the MHD numerical scheme) is essentially the same as that in 
Paper I and II. We assume an initial hydrostatic equilibrium in 
the gravitatio nally s tratified medium, adopting an iso-thermal 
King model (King 1962). The magnetic flux and the mass 
are steadily injected in a central small volume during a cer- 
tain time period. Since the injected magnetic fields are not 
force-free, they will evolve as a "magnetic tower" and inter- 
act with the ambient medium. 

2.1. Numerical Set Up 

In this paper, we present mainly two different cases: one 
is called the "unperturbed case" in the following discussion, 
which is a run without any initial perturbation to the back- 
ground initial profiles. This simulation is taken from the 
run presented in Paper II. The other is called the "perturbed 
case", where a finite amplitude perturbation (a few percent of 
the background sound speed) is given to the velocities of the 
background gas. 

The total computational domain is taken to be 
\x\, \y\: \z\ < 16. The numbers of grid points in the 
simulations reported here are 240 3 in the unperturbed case 
and 320 3 in the perturbed case, where the grid points are as- 
signed uniformly in the x, y, and z directions. We normalize 
physical quantities with the unit length scale Rq, the unit den- 
sity po, the sound speed C s q as the typical velocity field in the 
system, and other quantities derived from their combinations, 
e.g., the typical time to as Rq/C s o, etc. Normalizing factors 
are R = 5 kpc, C s0 = 4.6 x 10 7 cm/s, t = 1.0 x 10 7 
yrs, po = 5.0 x 10~ 27 g/cm 3 , the unit pressure po as 
PoCgQ = 1.4 x 10~ n dyn/cm 2 , and the unit magnetic field 

Bo as y/ AnpoCso = 17.1 /xG. We take po andpo as the initial 
quantities at the origin (x, y, z) — (0, 0, 0) and the p and 
p at the origin in the simulations are set to unity. The initial 
sound speed in the simulation is constant, C s = 7 1 / 2 w 1.29, 
throughout the computational domain. The simulation 
domain is from —80 to 80 kpc, with Ax = Ay = Az ~ 0.13 
for the unperturbed case and Ax = Ay = Az ~ 0.1 for 
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the perturbed case, which correspond to ~ 0.67 kpc and 
~ 0.5 kpc respectively. An important time scale in the 
system is the sound crossing time t s « 0.78, corresponding 
to a typical time scale Ro/C a o ~ 10.0 Myrs. Therefore, 
t = 1 is equivalent to the unit time scale 12.8 Myrs. 

In the King model we use here, we adopt the cluster core 
radius R c to be 4.0 (i.e., 20 kpc) and the slope k to be 1.0 in 
the unperturbed case and 0.75 in the perturbed case. The in- 
jections of magnetic flux, mass and its associated energies are 
the same as those described in Paper I. The ratio between the 
toroidal to poloidal fluxes of the injected fields is character- 
ized by a parameter a — 15, which corresponds to ~ 6 times 
more toroidal flux than poloidal flux. The magnetic field in- 
jection rate is described by 7& and is set to be 7& = 3. The 
mass is injected at a rate of j p = 0.1 over a central volume 
with a characteristic radius r p — 0.5. Magnetic fluxes and 
mass are continuously injected for tj n j = 3.1, after which the 
injection is turned off. These parameters correspond to a mag- 
netic energy injection rate of ~ 10 43 ergs/s, a mass injection 
rate of ~ 0.046Af Q /yr, and an injection time ~ 40 Myrs. We 
use the outflow boundary conditions at all outer boundaries. 
Note that for most of the simulation duration, the waves and 
magnetic fields stay within the simulation box, and all mag- 
netic fields are self-sustained by their internal currents. 

2.2. Definition of Forward and Return Currents 

As examined in Paper II, the magnetic tower jets can prop- 
agate without any visible distortion throughout the time evo- 
lution (see Fig. 1 of Paper II). The magnetic tower has a well- 
ordered helical magnetic configuration; a tightly wound cen- 
tral helix is going up inside the tower jet and a loosely wound 
helix is coming back at the outer edge of magnetic tower jet 
(see Fig. 3 of Paper II). The axial current flow J z displays 
a closed circulating current system in which one path flows 
along the central axis (the "forward" current) and another con- 
ically shaped path that flows outside (the "return" current) 
(see Fig. 4 of Paper II). Here, we define the forward current 
density J F and the return current density J R as follows; 



J R 



(Jz 
(Jz 



>0) 
<0) 



(1) 



Note that J F (and J R ) is a vector and can have additional 
radial and azimuthal components J F and J F (J R and J R ) 
with any sign. 

2.3. The Kruskal-Shafranov Criterion 

Accordin g to the well-known Kruskal-Shafranov (K-S) cri- 
terion ( Shafranov 1 119571 iRruskal et al. lfl958l) . a magnetized 
flux tube (along the z axis) is stable to the kink (m = 1) mode 
as long as the magnetic twist angle $(r) is below some criti- 
cal value $ cr it, 
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crit ; 
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where, L is the length of the current-carrying magnetic flux 
system, r is the cylindrical radius of the flux tube, and B z and 
Bcf, are the axial and azimuthal field components in a cylin- 
drical coordinate, respectively. In the original K-S criterion, 
$ cr it is equal to 2n. The effect of "line-tying", i.e., fixing the 
radial motion of the magnetic fluxes at the boundary, how- 
ever, could raise the stability threshold. (In the solar coronal 
loops, the foot points of magnetic fluxes can be anchored in a 
high density region as chromosphere.) The modified critical 



values for a force-free magnetic loop configurat ion in solar 
coronal loops are between 2tt and 107T (e.g.. lHood & Priest! 
ll979t lEinaudi & Van Hoven II1983T) . For the configuration of 
an astrophysical jet, which has L ^ r, inequality (0 can be 
easily violated. If the MHD jets originate from an accretion 
disk, the bottom part of jets may be line-tied in the disk and 
the top part of jets can be free against the radial motion; i.e., a 
"semi" line-tying configuration. (Further argument of the sta- 
bilizing effect by semi line-tying will be presented in §4.2.) 
For lack of detailed linear analysis of "semi" line-tying mag- 
netic configurations, we still define $ cl -it = 2tt, which should 
only be regarded as a guide for gauging the stability, rather 
than a precise threshold. Here, we replace L by a spectrum of 
axial wavelengths A, and define the critical wavelength to be 
the one that is marginally stable: 
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The K-S stability criterion (eq. |2)) then becomes 

A < A cr it, 



(3) 



(4) 



so that wavelengths greater than A cr ;t could be unstable to the 
CDI kink mode. 

2.4. Diagnostics of Non-axisymmetric Mode 

To better diagnose the non-axisymmetric distortion, we per- 
form the power spectrum analysis of the magnitude of the cur- 
rent density \J\. We analyze the modal structure of \J\ by 
performing the volume-averaged Fourier transform: 



/(to, k) 



I 

Vol 



\J\e i( - m4 ' +kz) rdrd4>dz, (5) 



■ 1/2 

where \J\ = (J? + Jj + J z ) and V c \ is the cylindrical 
volume that encloses J(= J F + J R ) 

fr b p2ir r z b 



r dr d(j) dz. 



(6) 



We use the values r a = 0.0, r\> = 7.75, z a = 0.0, and Zb = 
15.5 in the "upper" volume and use z a = —15.5 and Zb = 0.0 
in the "lower" volume. The quantity /(to, k) is a function of 
the azimuthal mode number to and axial wave number k — 
2ir/\ (corresponding to a characteristic wave length A), and 
it is a function of time as well. 
Finally, we identify the power spectrum as 

\f(m, k)\ 2 = {Re[/(m, k)]} 2 + {lm[/(m, fc)]} 2 ,(7) 



Re [/(to, k)] 
1 



VcU J J,, 
Im[/(m, k)] 
1 



| J| cos (to0 + kz) rdrd(j)dz, (8) 



| J\ sin (to0 + kz) r dr d(f>dz. (9) 



We will examine the time-dependent behavior of the power 
spectrum and the growth of each mode. 
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3. RESULTS 

3.1. Unperturbed Case — Dynamical Stabilization against 
the CDI Kink Mode 

We first re-examine the numerical result of Paper II as an 
unperturbed case, from a viewpoint of stability. For that simu- 
lation, we did not impose any initial perturbations to the back- 
ground state. So the only available perturbations are from the 
numerical noise in the double precision computations. Fig- 
ure ^ shows the snapshots of J z (top) and A cr j t (bottom) in 
the final stage (t = 10.0). Again, we can identify the cir- 
culating path of axial current flow J z which consists of both 
the forward axial current density J z ( J z > 0) and the return 
current (J z < 0) in the top panel ( jj lies around the cen- 
tral axis, while jf surrounds J z ). The maximum strength of 
| J~\ is much larger than that of | Jf~\, indicating that a tightly 
wound central helix is going up at the center of the tower and 
a loosely wound helix is coming back at the tower edge. The 
paths of forward and return currents are approaching one an- 
other around the mid-plane \z\ < 4.0. On the other hand, the 
path of return current is radially expanding at the upper axial 
position z > 4.0. This causes a separation of current flow 
paths in a thermally confined MHD jet with a helical mag- 
netic field. (This plays an important role in the excitation of 
the external kink mode and it will be discussed in §4.3.) 

As seen in the bottom panel, A cr it has a range 2.0 ~ 6.0 
in the core part of J z , and has no axial dependency except 
at the tower front. While A cr j t has no range in any part of 
meaning that the CDI kink mode is unlikely to grow at 
the return current region. Taking the whole part in the axial 
direction into consideration, the maximum length of current- 
carrying column is L ~ 16 at t = 10.0, which is larger than 
the distributed A cr jt in the core part of Jj. However, the non- 
axisymmetric displacement does not seem to occur as seen 
in both panels. We next show the result of the power spec- 
trum analysis of J F in order to confirm the non-existence of 
the growing kink mode. Figure|2]shows the space-time (A, t) 
diagrams of the power spectrum \ f(m, A)| 2 for m = 1 in 
the unperturbed case. We can not detect any growing non- 
axisymmetric mode (m = 1) over the wider range of axial 
wave length (A = 0.26 ~ 10.42). The time variation of 
the power spectrum of f(m, A) with a selected wavelength 
A = 6.0 in Fig. [3] Throughout the time evolution, the power 
in m = 1 mode is always below the axisymmetric (m = 0) 
one. Note that the variation of the m = mode (seen in both 
volumes around t = 4 ~ 8) is due to the axisymmetric defor- 
mation of the magnetic tower into a slender-shaped structure 
by the gravitational contraction radially. 

It may seem puzzling that, for the unperturbed case, the 
magnetic tower jet is stable even with L > A cr it. Here, we 
consider a possible mechanism of the stabilizing effect in our 
results. Figure 0] shows the time evolution of the K-S criti- 
cal wavelength A cl jt in the unperturbed case. A gradual in- 
crease of the distributed A cr it( = %TrrB z / B$) in the forward 
axial current-carrying column J z can be observed. This cor- 
responds directly to a decrease of the ratio B^/ B z . Because 
the injection of magnetic fluxes into the central volume in the 
system has been turned off at t = 3.1, the twisted magnetic 
fluxes are naturally relaxed during the expansion of the mag- 
netic tower. The threshold length of the marginally stable 
of J z in the axial (z) direction becomes longer as the mag- 
netic tower proceeds. Therefore this "dynamical relaxation" 
of magnetic twists plays an important role in the stabilization 
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FIG. 1. — Top: Distribution of the axial current density J z in the x — z 
plane at t = 10.0 in the unperturbed case (Paper II). Bottom: Distribution of 
the critical wavelengths A cr i t at t = 10.0 in the x — z plane, corresponding 
to the original Kruskal-Shafranov criterion for CDI m = 1 mode in the 
unperturbed case (Paper II). 



in our current results. We will discuss other effects on the 
stabilization in §4.2. 

3.2. Perturbed Case — Growth of the CDI External/Internal 
Kink Modes 

We next examine the non-axisymmetric evolution of the 
magnetic tower jet in a gravitationally stratified atmosphere. 
To excite non-axisymmetric modes (m > 1), a small random 
velocity perturbation (4% of initial sound speed) has been 
added into the system as an initial condition. 

Figure|5]shows a snapshot of the density p at the final stage 
t = 9.5 (see Fig. 1 of Paper II for comparison) in the per- 
turbed case. Even though the configuration of the magnetic 
tower is non-axisymmetric, the main features remain simi- 
lar to the unperturbed case. In general, the structure of the 
magnetic tower consists of both a well-collimated "body" and 
radially extended "lobes". Several key features in the mag- 
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FIG. 2. — Space-time (A, t) diagram of the azimuthal Fourier power (a natural log scale) for CDI mode m ■■ 
upper volume z > 0. Right: For the lower volume z < 0. 




1 in the unperturbed case (Paper II). Left: For the 



netic tower jets (the MHD wave structures, the heating pro- 
cess at the tower front, the cylindrical radial force equilibrium 
at the tower edge, and dynamically collimating process) are 
the same as the unperturbed case (see Paper II for details). 
The Alfven speed becomes large (about three times of local 
sound speed), while the plasma (3{= 2p/B 2 ) becomes small 
(j3 < 0.1) inside the density cavities due to the expansion of 
magnetic fluxes. The bulk flow inside the magnetic tower is 
sub-sonic and sub-Alfvenic, while the fronts of the preceding 
hydrodynamic shock and the magnetic tower can propagate 
super-sonically. 

Similar to Fig. |2] we present the space-time (A, t) diagrams 
of the power spectrum \f(m, A) | 2 in the perturbed case for the 
m = 1 and 2 in Fig. [6] Clearly, the growth of the m = 1 mode 
can be seen over a wide range of A after t ~ 4 in both upper 
and lower volumes. Furthermore, the m = 2 (elliptical) mode 
also grows at a relatively short wavelength around A = 2 ~ 3 
after t ~ 6. To identify the growth rates of different modes, 
we show the time variation of the power spectrum of f(m, A) 
in the perturbed case with a selected wavelength A = 6.0 in 
Fig. ^(similar to Fig. [3}. The m = 1 mode in both volumes 
exhibits an approximate exponential growth on a dynamical 
timescale during t = 4 ~ 6. At later stages, after the m = 1 
mode grows above the axisymmetric m — mode, its growth 
rate decreases slightly, but still keeps growing exponentially 
to saturation at the final stage. The m — 2 mode also shows 
the similar behavior with the m = 1 mode at the earlier stage 
around t = 6 ~ 8, but it saturates at or below the level of the 
m = mode. 

A linear fit using a least-square method has been performed 
on the initially exponential growth for the m = 1 mode to 
extract a linear growth rate 



Im (u>) 



dln\f(m,X)\ 
dt 



(10) 



The estimated growth rates in the perturbed case are also 
given in Fig. The growth rate in the upper volume ap- 
pears to be a bit stronger than that in the lower volume. We 
confirm the distributed Alfven speed Va in the surfaces of J F 
is of the order of 8 — 9 in our normalized unit during the time 
evolution. The spatial length scale £ of the non-axisymmetric 
(m > 0) distortions in the radial direction is of the order of 
2 — 4, as seen in Fig. |5] Therefore, the inverse of the Alfven 
crossing time r^ 1 in the distorted part of magnetic tower jet 



is given by 



Va 



2.0-4.5. 



(11) 



This is, in general, consistent with the timescales of the grow- 
ing m = 1 mode derived from Fig. Q So, our results lead us to 
conclude that the non-axisymmetric distortion in the magnetic 
tower jet is caused dominantly by the normal current-driven 
kink (m = 1) mode. 

Two types of CDI modes are identified in the context of 
magnetically cont rolled fusion pla sma; the "external" and the 
"internal" modes (Batema n 119781) . The external mode can be 
identified with a radial distortion, while in the internal mode, 
there is no radial perturbation at the outer edge of the current- 
carrying column. One of our findings in the perturbed mag- 
netic tower jet is that the CDI external and internal modes "co- 
exist" during the dynamical evolution. Figure [8] shows snap- 
shots of J z in the perturbed case at t = 8.0 and 9.5. As seen 
in the top panel (t = 8.0), the edge of J F is clearly shifted 
to an offset axial direction at higher latitude {\z\ > 4). On 
the other hand, that edge never shifts non-axisymmetrically 
at lower latitude (\z\ < 4). However, the non-axisymmetric 
distribution of J F at the inner part can be detected, indicat- 
ing that the internal mode is certainly growing. At the final 
stage (t — 9.5: bottom), we can see the well-grown wiggles 
inside the "lobes", while the edge of the "body" part in the 
magnetic tower still remains quasi-axisymmetric except for 
the inner part. Taking the result of Fig. |6]into consideration, 
the m — 1 kink mode is dominant in the external CDI at the 
lobes, while the m = 2 elliptical mode may be dominant in 
the internal CDI at the body of the magnetic tower. In addi- 
tion, we illustrate a 3-D view of some selected magnetic lines 
of force in the perturbed case at the final stage t = 9.5 in Fig. 
[9] It indicates that a tightly wound central helix, which goes 
up along the central axis, relaxes in the lobe part where the 
external kink mode has taken place, and a large-scale loosely 
wound twist is formed there. 

4. DISCUSSION 

4.1. Suppression of Shear-driven Instability 

We first discuss the possibility of the interplay between CDI 
and KHI modes in our results. The examinations of KHI in 
linear/nonlinear regimes have been extensivel y performed by 
P. E. Hardee and his collaborators (see, e. s. . lHardee 1120041 
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FIG. 3. — The time variation of the azimuthal Fourier power (a natural 
log scale) in modes m = 0, 1 with a specified wave length A = 6.0 in the 
unperturbed case (Paper II). Top: For the upper volume z > 0. Bottom: For 
the lower volume z < 0. 



and references therein). In general, non-axisymmetric "sur- 
face" (external) modes (m > 0) are stable against KHI for 
sub-Alfvenic flow. However, in super-Alfvenic but trans- 
fast magnetosonic flow, they can be unstable. Only ax- 
isymmetric surface mode (m = 0) remains unstable in sub- 
Alfvenic flow, alth ough with a relatively small growth rate 
iBodo et al. II1989I) . A potential zone of unstable region in 
non-axisymmetric KHI surface modes should exist in down- 
stream of the Alfven surface. 

The instability crite rion for non-axisymme tric KHI surface 
modes is as follows (Hardee & Rosen 2002) : 



PP.. 




> V As 



I Pj +Pe 



(12) 



FIG. 4. — Time evolution of the K-S critical wavelength A cr i t in the x — z 
plane at t = 5.0 (Top), 7.5 (Middle), and 10.0 (Bottom) in the unperturbed 
case (Paper II). 



where AV = \Vj — V e \ is the velocity shear and Va s is the 
surface Alfven speed (subscript j corresponds to the jet itself 



Non-axisymmetric Magnetic Tower Jets 



7 



a I—, 



— J.O 



V: 
3 



■■_ 

"J 




FIG. 5. — Distribution of the density p (logarithmic scale) in the x — z plane is shown along with the poloidal velocity field (arrows) at t = 9.5 in the perturbed 
case. Similar to the unperturbed case (Paper II), density cavities are formed due to the magnetic field expansion, while the non-axisymmetric motion can be seen. 
The length of the arrow at the upper right shows the unit speed. 



and e corresponds to the external medium). On the other hand, 
the non-axisymmetric "body" (internal) modes (m > 1) can 
become important and affect the jet interior in the following 
situations: (1) if the jet velocity exceeds the super-fast mag- 
netosonic speed Vf. 



Vj > v f , 



(13) 



or (2) if the flow velocity is slightly below the slow magne- 
tosonic speed V s : 



C S V A 



< Vj < v., 



(14) 



where C s is the sound speed and Va is the Alfven speed 
(Hardee & Rose n 119991) . 

Figure EH shows the transverse profiles in the x direction 
of several quantities at t = 8.0 in the perturbed case. Two 
different axial positions z ~ 2.0 and 8.0 are selected for ex- 
amining the KHI surface/body modes. For the surface modes, 
we inspect the transverse distribution of the bulk flow speed 
V, the density p, and the magnetic field strength B. Two dis- 
tinct velocity shears (in both positive and negative x regions) 
are identified in the x direction on z = 2.0 (upper-left); an 
outer part 0.6 < \x\ < 1.5 and an inner part 0.2 < \x\ < 0.6. 



At the outer part, AV ~ 1.2 and Va s ~ 1-8 and at the inner 
part, AV ~ 1.1 and Va s ~ 8.7. Therefore, the inequal- 
ity AV < Vas is satisfied. Similarly, we can see many dis- 
tinct velocity shears in the upper-right panel. However, even 
in these shears, we confirm that the inequality AV^ < Vab 
still holds everywhere. This indicates that the KHI surface 
modes are completely stabilized, because there is no signifi- 
cant velocity shear exceeding the local surface Alfven speed 
in the transverse (x) direction. For the body modes, we in- 
spect the transverse distribution of the bulk flow speed V, the 
Alfven speed, and the quantity C s Va/{C^ + Vl) 1/2 . Inside 
the magnetic tower, the plasma (3 (— 2p/B 2 ) is always less 
than unity during the time evolution. Therefore, the phase 
speed of the MHD fast mode wave is equal (in the local mag- 
netic field direction) or larger (in the oblique direction of the 
local magnetic field) than the Alfven speed, Vf > Va- As 
seen in both Bottom panels in Fig. ^0| the bulk flow speed is 
always extremely sub-Alfvenic, and furthermore, the inequal- 
ity V < C S V A /(C 2 + Vl) 1 ' 2 is also satisfied in the "body" 
part of the magnetic tower jet. These indicate that the KHI 
body modes are suppressed as well. We conclude that the 
KHI modes play no role and the normal CDI modes create the 
non-axisymmetric structures in our results. 
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Fig. 6. — Similar to Fig. |2] but for modes m = 1 {Top) and 2 (Bottom) in the 
z < 0. 

4.2. Stabilization/De stabilization of Magnetic Tower Jets 

As examined in A3. II the unperturbed magnetic tower jet 
can be stabilized against the CDI kink mode beyond the K- 
S criterion. One of the certain mechanisms of the stabilizing 
effect, which we can identfy in our unperturbed case, is the 
dynamical relaxation of magnetic twists. However, this pro- 
cess may contribute the marginal stabilization under an ideal 
situation (no perturbation except for the numerical noise is as- 
sumed as an initial condition). If the system is perturbed ran- 
domly with a finite amplitude (a few percent of the local sound 
speed), then the non-axisymmetric (m > 1) internal/external 
CDI modes can grow. The disturbances in the system are rea- 
sonable (especially when related to possible astrophysical sit- 
uations). Therefore, the magnetic tower jet in our treatments, 
will eventually be subject to CDI. Here, we give further dis- 
cussions on the stabilities of magnetic tower jets. 

4.2.1. Line-tying Effect 

In the solar coronal magnetic loops, it is well known that the 
line-tying effect at the foot points (chromosphere) can raise 
the original K-S stability threshold (dense materials can pre- 
vent the radial motion of twisted magnetic fluxes). How about 
the case in the astrophysical jets? Anchored magnetic fluxes 
to the accretion disk and/or stellar/black hole surfaces may be 
considered as a semi (one side) line-tying boundary condition 
in a certain dynamical stage. An important clue to this can 
be derived from laboratory experiments of magnetic towers 
JHsu & Bellan ll200l 120031 Their quasi "force-free" mag- 
netic towers can exhibit the systematic m = 1 distortion when 
the current-carrying system satisfies the original K-S limit. It 
might be unlikely that such a semi line-tying effect can stabi- 
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perturbed case. Left: For the upper volume z > 0. Right: For the lower volume 

lize the jet even in the astrophysical situation. In our global 
modeling of magnetic tower jets, we do not solve the evolu- 
tion of the central region (an accretion disk, etc.). Therefore, 
we can not apply it directly to the AGN central region. It is 
still an open issue whether a semi line-tying effect can work 
or not in the real astrophysical jets. 

4.2.2. Internal Thermal Pressure 

As discussed in Paper II, the "forward" current-carrying 
column J F never reaches a force-free equilibrium. Instead, 
it will be supported by the outward-directed, local (core of 
the magnetic towers) pressure gradient force: — Vp + J x 
B ~ during a dynamical evolution (see Fig. 10 in Pa- 
per I I). In terms of m inimization of the MHD energy inte- 
gral iBernstein et al. lfl958l) . a finite gas pressure in the sys- 
tem (plasma (3 ^ 0) can contribute to the stabilization. How- 
ever this can be effective when /3 > 1, i.e., the gas pressure 
is more domina nt than the magnetic pressure (for m = 1 
kink mode, see lHood&Priestlll979t) . Therefore, we may 
not expect a strong stabilization by the pressure in the mag- 
netic tower ((3 < 0.1). Furthermore, the outward-directed 
gas pressure gradient itself (— Vp > 0) actually destabilizes 
the tower through the PDI (interchange and/or ballooning) 
modes, which is driven by the perpendicular currents (in con- 
trast to the axial [parallel] currents, considered mainly in this 
paper). We, however, could not find any growth of the PDI 
modes, in both the unperturbed (Paper II) and the perturbed 
cases, implying no direct contribution of the pressure gradient 
to the destabilization by the PDI modes (see next section for 
more details). Also, we can not confirm any direct contribu- 
tions of the internal thermal pressure effects (fi ^ and — Vp) 
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Fig. 7. — Similar to Fig.l3l but for modes m = — 2 in the perturbed case. 
Straight lines are derived by fitting the slopes of the m = 1 mode growth to 
a linear function. Top: For the upper volume z > 0. Bottom: For the lower 
volume z < 0. 

to the destabilization by the CDI internal/external modes. 

4.2.3. Magnetic Field Configuration 

We next consider the magnetic field configuration in mag- 
netic tower jets. In the context of magnetic controlled fusion 
systems, the helical field configuration in the tower model can 
be identified as the reversed field pinch (RFP) profile as seen 
in Fig. 1 1 in Paper II. The axial component B z has a crucial 
shear (its sign reverses from "+" to "-") in the radial direction 
(RFP in the context of astrophysical jets, see also Benford 
2006). Several stab ilization effects in RFP are summarized 
in FreidbgraHilll; (i) a hollow or very flat pressure profile 
(Vp ~ 0) are required near the central (z) axis to suppress the 
PDI sausage (interchange, internal) mode (to = 0), but (ii) it 
can be suppressed when (3 < 0.5. This may be the situation 
in the presen t paper. As already discussed in the previous sec- 
tion ( §4.2.2\ . a finite pressure gradient exists to contribute to 
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FIG. 8. — Snapshots of the axial current density J z in the x — z plane at 
t = 8.0 (Top) and 9.5 (Bottom) in the perturbed case. Two types of growing 
kink modes ("ExternarV'Tnternal") are illustrated on the Top panel. 

the radial force equilibrium in the core part (around the cen- 
tral axis) and this might play a role in driving the PDI modes. 
However, the plasma (3 is small (/3 < 0.1), indicating that the 
destabilization by the PDI modes may be suppressed due to 
the small plamsa j3. Also, (iii) the strong B z shear can prevent 
the destabilization by the CDI internal/external kink modes 
(to = 1). In our results, this shear may provide the stabiliza- 
tion against the CDI external modes, while it does not affect 
the CDI internal modes at the lower latitude \z\ < 4, as seen 
in Figs. |6]and|8] Furthermore, the CDI external modes can 
grow at the higher latitude \z\ > 4. This may be because the 
strong B z shear in the radial direction will be relaxed towards 
the jet axis due to the decreasing external pressure (expanding 
the poloidal fluxes to the radial direction). 

4.2.4. External Environment 

The external gas itself also may affect the stabilization. 
The inward-directed, thermal pressure confinement and/or 
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FIG. 9. — Three-dimensional view of some selected magnetic lines at 
t = 9.5 in the perturbed case. 



the gravitational contraction of the external unmagnetized 
gas, which surrounds the magnetic tower jets, may provide 
some stabilization at least against the growth of external CDI 
modes, although the internal CDI and PDI modes will not 
be affected (for the i nternal PDI sausage m = mode, see 
Lerv & Frank 2000). In our results, the external environ- 
ment may have a similar stabilizing effect on the external 
CDI modes to the B z shear, as discussed in the previous sec- 
tion ( fl4.2.3> . while the internal PDI mode may be suppressed 
by the low internal pressure. In the magnetic tower jets, an 
external gas plays a crucial role in the stabilization of the 
CDI external modes. When an external gas pressure profile 
does not decrease (constant), which is examined in Paper I 
(see Fig. 1 1 in Paper I for a continuous injection case), the 
magnetic tower never expands radially and both the thermal 
confinement by the external gas and the strong B z shear will 
be kept throughout the magnetic tower "body" (no formation 
of "lobe"). Therefore, the CDI internal modes may occur, 



but the CDI external modes should be surpressed . These are 
consistent with the CDI "non-disrup ting" model (lApplll996t 
Appl et al. l2000tlLerv et al. 1 2000) as well. 

4.2.5. Rotation 

We finally discuss the stabilization by the rotation. In our 
numerical treatment, we investigate almost "non-spinning" 
magnetic tower jet. From the theoretical point of view, 
any magnetically driven outflows (magneto-centrifugal and/or 
magnetic pressure driven) can have an azimuthal velocity 
component. They possess an angular momentum provided 
by the accretion disk or Kerr black hole magnetosphere. In 
the classical model of nonrelativistic MHD winds, the cen- 
trifugal force plays a minor role in the radial force e quilib- 
rium far from the accretion disk (Blandford & Pavne 1982). 
However, this is only the case in the nonrelativistic treatment. 
In the relativistic MHD regime, the situation will be modi- 
fied; semi-analytical solutions of relativistic MHD winds with 
self-similarity exhibit that the relativistic inertia and electro- 
static forces play an important role in the radial force equilib- 
rium at the cylindrical flo w struct ure, extended to pc or larger 
scales ( Vlahakis & Konigl 2003). Braided emission-line pro- 
files and rotational velocity in NGC 4258 indicate a pure heli- 
cal mo tion exists alon g the tightly wrapped strand even in kpc 
scales (ICecil. Wilson. & Tullv 11991 . 

The stabilizing effect of rotation against the CDI kink mode 
in the spinning nonrelativistic PFD/KF D jets has been inves - 
tigated by numerical simulations (Nakamura & Meier 2004). 
Because of the centrifugal effect, these spinning jets can be 
stabilized against the m = 1 mode beyond the point pre- 
dicted by the K-S criterion. The linear stability analysis 
of the relativistically rotating force -free fields in Kerr ge- 
ometry (Blandford & Znaiek 1977) has been performed by 
Tomimatsu. Matsuoka. & Takahashi (2001). They also show 
that the field-line rotation has a stabilizing effect against the 
kink modes satisfying the K-S criterion. In these cases, the 
presence of an azimuthal velocity component can modify the 
radial force equilibrium, making the jet robust against the CDI 
kink mode. Examinations of the stability property in a rotat- 
ing magnetic tower jet are clearly desirable. 

4.3. External CDI Kink Mode as a Possible Mechanism for 
Wiggling AGN Jets 

The CDI of current-carrying jets for two types of re- 
turn current distributions have been discussed: (1) the 
return current is assumed to flow within the jet itself 
(the jets are thermally confined by the extern al medium) 
dChiuderi. Pietrini. & Torricelli-Ciamponi 119 89). and (2) the 
return current is assu med to flow around t he jet, such as in a 
magnetized cocoon (Benford 1978, 2006). Most of the theo- 
retical/numerical works have been devoted to the former case, 
in which the forward current J F and the return current J R are 
concentrated at the interfaces between the jet and ISM/IGM, 
which is taken to be "unmagnetized". In this situation, the 
growth of external CDI modes can be effectively suppressed, 
while the excitations of only internal CDI modes are allowed. 
There is no radial displacement at the edge of the current- 
carrying column by the internal CDI modes, and therefore, it 
may be difficult to explain the off-set axial distortion of ob- 
served jets by these internal modes. As described in fJO this 
is one of the reasons why much less attention has been paid 
to the CDI than to the KHI for the stability problems in astro- 
physical jets to explain the disrupted jet exterior. 
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FIG. 10. — Transverse profiles in the x direction of several quantities at t 
The bulk speed V ( 



8.0 on 2 = 2.0 (Left panels) and on z = 8.0 (Right panels) in the perturbed case. 



•J V£ + Vy + V§), the density p, and the magnetic field strength B (= B% + + Bf) are shown for inspecting the KHI surface 

modes (Top panels). The bulk speed V, the Alfven speed Va, and the quantity C s Va/(CJ + V^) 1 / 2 are shown for inspecting the KHI body modes (Bottom 
panels). 

However, in the gravitationally stratified atmosphere, the 5. CONCLUSION 

effect of thermal confinement on the edge of the current- 
carrying column will get gradually weaker as the magnetic 
tower grows. As a result, a separation of current flowing paths 
between J F and J R occurs, as seen in Fig. [8] This implies 
that the edge of J F becomes a free boundary against the CDI 
external modes, even when it is embedded deeply inside the 
external, unmagnetized gas. Furthermore, the region between 
the outside of J F and the inside of J R becomes nearly force- 
free J x B ~ 0. Therefore, the CDI external modes can 
grow in a dynamical time scale to form the systematic wig- 
gles. This is also true even in the spinning PFD/KF D jets 
under the large-scale mag netic field ( Nakamura et al. ll200lt 
Nakamura & Meier 2004). On the other hand, even if the gas 
pressure in the external, unmagnetized gas is high enough to 
suppress the growth of the CDI external kink mode, then the 
CDI internal kink mode can grow to form the internal non- 
axisymmetric structure, as we exhibited in Fig. [8] Therefore, 
we believe that the external CDI modes are one of the possi- 
ble mechanisms to explain the disruption of current-carrying 
MHD jets and it could be applied to the wiggling AGN jets. 



Stability properties of magnetic tower jets have been exam- 
ined by performing 3-D MHD simulations. Our numerical 
results show that the magnetic tower jets in the gravitationally 
stratified atmosphere could survive the current-driven kink in- 
stability beyond the Rruskal-Shafranov criterion even in the 
nonlinear regime. One of the potentially stabilizing effects 
on a traditionally unstable magnetic field may be due to the 
dynamical relaxation of magnetic twists in the propagating 
magnetic tower jets. This causes an increased threshold of the 
unstable critical wavelength. 

However, the magnetic tower jets are eventually distorted 
by the non-axisymmetric perturbations with a few percent of 
the local sound speed in the system. The current-driven kink 
modes grow predominantly on time scales of the order of the 
local Alfven crossing time ta- In the gravitationally stratified 
atmosphere, two types of kink modes appears: the internal 
and external modes. At a large distance away from the central 
region, the external kink mode grows, while only the internal 
kink mode exhibits near the central region. Therefore, the 
exterior of magnetic tower jets will be deformed into a large- 
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scale wiggled structure by the external kink mode. 

Possible mechanisms of stabilizing/destabilizing in the 
magnetic tower jets under the gravitationally stratified atmo- 
spheres are discussed. None of the growing surface/body 
modes of Kelvin-Helmholtz instability are identified due to 
the small velocity shear and slow bulk speed (compared with 
the local Alfven speed). The pressure-driven instabilities are 
also inhibited due to the low plasma (< 0.1) inside the 
magnetic towers. The "reversed field pinch" field configu- 
ration in the towers and the external unmagnetized gas may 
provide some stabilization effects on the current-driven exter- 
nal modes. 

Non-axisymmetric current-driven instabilities (m > 1) 
are absolute instabilities, i.e., they grow but do not prop- 
agate in the jet co-movin g (rest) frame (IAppI etal.lbOQOt 
Nakamura & Meier 2004). If the dissipation timescale t& of 
the current associated with jets is much longer than the jet 



dynamical timescale Tj < ta ~ 7"cdi <C t^, the patterns cre- 
ated by the CDI modes will persist for some time as the flow 
advances. And the bulk flow itself will appear to travel in a 
true 3-D helical pattern as it follows the magnetic backbone of 
the helix, as obse rvationally expected in several radio source s 
as 3C 345 JZensus etal. 11 995). 3C 120 JGomez etal. 120011) . 
and 3C 449 (F eretti et al. 119991) . 
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